Critical Phenomena: An Introduction from a modern perspective * 
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Our aim in this set of lectures is to give an introduction to critical phenomena that emphasizes the emergence of and 
the role played by diverging length-scales. It is now accepted that renormalization group gives the basic understanding 
of these phenomena and so, instead of following the traditional historical trail, we try to develop the subject in a way 
that emphasizes the length-scale based approach. 
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I. PREAMBLE 



A phase transition is defined as a singularity in the free-energy or any thermodynamic property of a system. "A 
system" is to be understood in a general sense; it is defined by a Hamiltonian under various external macroscopic 
constraints. For example, an isolated collection of particles would be an example of a system, just as the same collection 
exchanging heat or work with surrounding. It will also be the same system even when it exchanges particles with 
surroundings. The three examples given are the standard micro-canonical, canonical and grand canonical ensembles 
in statistical mechanics. There are many other types of ensembles. All these ensembles are expected to give the 
same description of the system in the thermodynamic limit, at least that's what says the principle of equivalence of 
ensembles. 

Phase transitions are generally classified as first order, second order ...etc, or first order and continuous type. The 
behaviour of any system near or at a continuous transition (also called critical point) is distinctly different from the 
behaviour far away from it. This peculiarity near a critical point is known as critical phenomenon. 

Power laws are distinctive features of critical phenomena and can be ascribed to diverging length-scales. This 
connection is now so well-characterized that any phenomena, equilibrium or nonequilibrium, thermal or nonther- 
mal, showing power laws tend to be interpreted in the same fashion of equilibrium critical phenomena through the 
identification of relevant terms, scaling and diverging lengths. 

In these notes, we shall concentrate on the general aspects of critical phenomena. A few features of the first order 
transition will be touched upon at some point. 

One of the most important contributions of the studies of critical phenomena is the shift in point of view to a 
length-scale based analysis and classifying terms of a Hamiltonian (or in any other description of a problem) as 
relevant or irrelevant rather than classifying them as numerically weak or strong. The approach we take focuses on 
this and related issues. 

Our approach is to show that the class of singular behaviour seen near criticality can be characterized and understood 
via an emergent diverging length scale as the transition point is approached. This way of characterizing singularities 
find straight-forward generalization to many phenomena beyond equilibrium phase transitions. 

Defining via singularity is the most general way of characterizing a phase transition. For a large class of systems, 
singularities could occur due to "ordering" after a phase transition (or symmetry breaking) but that is not necessarily a 
requirement for all transitions. In other words, the existence of an "order-parameter" like quantity is not a prerequisite 
for understanding criticality though it might be useful in many contexts. We therefore deviate from the conventional 
historical approach to the subject via mean field theory, Landau theory, experimental observations leading to scaling. 

Outline of the notes: The problem from a theoretical point of view is presented, with a little bit of recapitulation 
of familiar things, in Sec II to IV. A possible resolution through scaling is introduced that leads to finite-size scaling. 
The consequences and physical interpretation leading to length-scale dependent parameters are discussed in Sec. V 
and VI. The Gaussian model and the 4> 4 theory are used as examples in section VII. The problems (exercises), even if 
not attempted, should be read for continuity. Only elementary calculus is used throughout, but the reader is assumed 
to have the background of statistical mechanics at the level of Reif. 

A. Large system: Thermodynamic limit 

In the canonical ensemble, the properties of the system with Hamiltonian H are obtained from the free energy, 

F = -k B T In Z, (1) 

where ks is the Boltzmann constant, T the temperature and Z = ^ e^ l3H is the partition function with (3 = 1/ksT. 
If a phase transition is a singularity of this free energy, then, for a simple well-defined H, it can occur only if Z has 
zeros for real values of the parameters. This cannot be the case for a finite-sized system because Z is a sum of a 
finite number of Boltzmann factors (which are always positive) . The zeros of Z or singularities of F can only be for 
complex values of the parameters. These complex zeros might have real limits when N — > oo ( infinite number of 
terms to be added in Z) . Then and only then can a phase transition occur. 

Conclusion: A phase transition cannot take place in a finite system. The thermodynamic limit, N —* oo, 
has to be taken. 
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II. WHERE IS THE PROBLEM? 



Statistical mechanics starts with the idea of a micro-canonical ensemble: the Boltzmann formula S = ks hif2 where 
O is the degeneracy or the number of states accessible to the system under the given conditions of energy (E) , volume 
(V), number of particles (N), etc. 

Thermodynamics starts with the basic postulate of the existence of an 'entropy' function S, so that 

E = E(V,S,N), or S = S(E,V,N). (2) 

One may consider either a system of fixed number of particles in equilibrium with surroundings or a large isolated 
system and focus on a small part which is in equilibrium with the rest. By releasing constraints (like E =const, or 
V = const, or N = const etc) one changes the ensemble. This is tantamount to changing the thermodynamic potential 
required to describe the system. 

Mathematically, the change of one thermodynamic potential to another is achieved via Legendre transformation. 
For example, from entropy to temperature (micro-canonical to canonical) is done via 

T= dE(V,S,N) 



( 3a ) 

F = E-TS, (3b) 

where F is now a function of (V, T, N). This is possible if and only if Eq. ([5a]) can be inverted to eliminate S on the 
rhs of Eq. |3^ in favour of T. This inversion is possible if dT/dS ^ 0. Similarly, we require dp/dV ^ 0, dh/dM ^ 0, 
or dpLj 'ON ^ etc, where /x is the chemical potential, h magnetic field and M magnetization (for a magnetic system), 
for change of appropriate ensembles. 

Since ^ is related to the heat capacity of the system, we see trouble if the specific heat diverges for some values 
of the parameters of the system. This looks like an algebraic problem of the transformation but it is vividly reflected 
in the argument for the equivalence of ensembles. So let us recollect that. 

Take the case of canonical and micro-canonical ensembles. The two ensembles are equivalent if the energy fluctuation 
in the canonical case is very small. First note that, if F oc N, (E) and (E 2 ) - (E) 2 are both proportional to N (more 
on this later). The angular brackets indicate statistical mechanical averaging. A straight forward manipulation, using 
derivatives of the partition function, yields 



where AE is the standard deviation of E and cy is the specific heat (heat capacity per particle). For a thermodynamic 
system, N — * oo, and this fluctuation goes to zero, provided c v is finite. Hence the equivalence. 
Well, the argument fails when cy — > oo. 

By our definition of phase transitions, a diverging second derivative of the free energy implies a phase transition. 
A hypothetical possibility? No, it occurs in many systems in many different contexts. Such points will be classified 
as critical points. Wait for a better definition of a critical point. 

Critical points seem to be at odds with the conventional wisdom of thermodynamics and statistical mechanics. 
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III. RECAPITULATION - A FEW FORMAL STUFF 



A consistent thermodynamic description requires a few postulates, some important ones of which are: 

(i) Existence of "entropy" ( already mentioned) 

(ii) Extensivity of E and S. 

(iii) Convexity of the free energy. 

Postulate (ii) is a simple additivity property, while (iii) is required for stability. The formulation of statistical mechanics 
guarantees (iii) but not (ii). We need to ensure (ii) to make contact with thermodynamics and this restricts the form 
of the Hamiltonians we need to consider. (See exercises III. 6 and III. 7 below). 



A. Extensivity 

From our experience of macroscopic system, we expect that under a rescaling V — > bV, N — > bN, and S — > bS, 
total energy should change accordingly, i.e. 

E(bV, 65, bN) = bE{V, S, N) (homogeneous function). (5) 

By choosing b — 1/N, (b is arbitrary), we have 

E{V,S,N) = N £(V/N,S/N), where £ (x, y) = E(x, y, 1). (6) 

Note that £ (v, s) is the energy per particle and being a function of v, s ( volume and entropy per particle), £(v, s) 
is independent of the overall size (e.g. N) of the system. This proportionality to the number of particles (or volume) 
is extensivity. 

Exercise III.l Additivity means E(Vi, Si, Ni) + E(V 2 , S 2 , N 2 ) = E{V 1 + V 2 ,S 1 + S 2 ,N 1 + N 2 ). Show that 
this implies Eq. [J. 

Exercise III. 2 Sanctity of extensivity: Remember Gibbs paradox and its resolution? 
The classical thermodynamics is based on Eq. || which, via the Euler relation for homogeneous functions, gives 

4f + "f + 4f = E = rs-™, (7) 

defining T = dE/dS, —p = dE/dV and fj, = dE/dN, where the partial derivatives are taken keeping the other 
variables constant. Each of the three new parameters defined, T,p and fi, is a derivative of an extensive quantity with 
respect to another extensive variable. Therefore T,p,fi, and variables like these are independent of the scale factor b 
of Eq.(3). Such quantities are called intensive quantities. Q 

Exercise III. 3 Take b = 1 + 81 and derive Eq. \%for 81 — > 0. 

Exercise III. 4 From Eq. M) show that 

F(bV, T, bN) = bF{V, T, N), and F(V, T,N) = N f(v, T). (8) 



1 Quantities like v, or s, an extensive quantity per unit volume or per particle, are also independent of the scale factor 6. To 
distinguish these from the other type, we may call them "densities" or "fields" (highly confusing for field theorists). In most 
cases, these could be distinguished from the context. 



4 



B. Convexity: Stability 

We recognize that the derivatives needed for the Legendre transforms (see Eqs. ([3a])) are derivatives of conjugate 
pairs: specific heat, compressibility, susceptibility etc. These derivatives (second derivatives of free energy) are called 
response functions, because they measure the change in the extensive variable as the externally imposed conjugate 
intensive quantity is changed. 

A generalization of the derivation of Eq. ^ for a pair (y, $) ($ being the extensive variable and y the conjugate 
intensive variable) leads to 

(9) 

which connects the response function for $ to the latter's variance. The latter being a positive definite quantity 
ensures that the response function is of a particular sign only. This in turn shows that the free energy as a function of 
y can have only one type of curvature (positive or negative) - a property known as convexity of the free energy. This 
positivity of the response function guarantees thermodynamic stability. This is the third postulate of thermodynamics 
mentioned above (quite often stated as the maximization/minimization principle). This connection between response 
and fluctuation plays a crucial role in the subsequent development, especially in developing a correlation function 
based approach. 

The main points we need here are 

(1) thermodynamic potentials are additive and therefore obey a simple scaling. 

(2) There are conjugate pairs of variables (T,S), (p,v), (/i, N), etc where the first one in each set is intensive (i.e. 
independent of the size) while the second one is extensive (i.e. proportional to the size). One may change variable 
(e.g., by releasing a constraint) and this is the change of ensemble in statistical mechanics. Two bodies in contact in 
equilibrium need to have equality of the relevant intensive quantities (remember zeroth law?). 

1. Comments 



are responsible for the problem with naive exten- 
sivity. 

• Let us write $ = j d d x </>(x), where </f>(x) is a lo- 
cal quantity (density). We have written it as a 
d-dimensional integral, though it would be a sum 
for discrete systems. The width of the proba- 
bility distribution for <f> around the mean (<£>) is 
given by the corresponding " susceptibility" as given 
by Eq. (|^). Assuming translational invariance, 
(<p(x)) = {4>), and denoting the response function 
per particle by %■> the relative width of P($) is given 
by N-Wy/xH®. 

Exercise III. 5 A very important model that is used extensively is the Ising model. 

H = -J^2s i s j -h^2s i , J,h>0, (10) 

<ij> i 

where the spins Si = ±1 are situated on a d-dimensional lattice, the interaction could be restricted to 
nearest-neighbours only (denoted by < ij >), and h is an external magnetic field. 

(1) Show that the free energy is extensive for any value of T and h, except possibly a particular point. 
(Difficult). 

Show this extensivity at T — and T — oo. (easy) 

(2) How do we define the dimensionality of the lattice? (Hint: How does the number of paths connecting 
two points grow with their separation?) For uniform hyper-cubic lattices (linear, square, cubic,...) d can 
be uniquely defined by the number of nearest-neighbours. 



dy 



k B T 



(($ 2 



• Equivalence of ensembles relies on the sharpness of 
the probability distribution for say E. The width of 
the distribution is related to the corresponding re- 
sponse function. For sharp distributions, the prob- 
ability of E for the combined system can be taken 
as the product of the individual probabilities (i.e. 
P Vl+ v 2 {Ei + E 2 ) = P Vl {E 1 )P V2 (E 2 ), whence fol- 
lows extensivity. Broad distributions would create 
problem here. 

• Broad distributions imply large fluctuations. It is 
therefore expected that fluctuations would play an 
important role in critical phenomena. Fluctuations 
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Exercise III. 6 Consider now more general Hamiltonians where each spin interacts with every-other: 
-Hbogus = — s i s j> an d #mf = — 5-/ S ' Sj ' ^ 1S tota l num ber of spins.) 

Show that there is a problem with extensivity for -ffbogus but not for H m f (for N — > oo^. Do this for T = 0. 
This model (H m f) gives the mean-field theory as an infinite dimensional model, -ffbogus is to be dumped. 



Exercise III. 7 Diamagnetism shows negative susceptibility. Any problem with convexity? 



IV. CONSEQUENCES OF DIVERGENCE - PROBLEM WITH EXTENSIVITY? 



Choose a quantity say the specific heat or susceptibility, that diverges at the critical point T — T c in the thermody- 
namic limit. For simplicity^] we keep only T as the control parameter (all others being kept at their respective critical 
values). Let us denote this quantity by C(T,N) - a total quantity for the whole system - obtained from the partition 
function of an TV-particle system. 

Extensivity requires that N~ 1 C(T, N) has an iV-independent limit for N — > oo, say c(T), so that for large but finite 

N 



C(T,N) = N c(T) + C cm (T,N), 



(12) 



Here C COI (T,N) is the correction term in an asymptotic analysis. An obvious expectation is ^■C cor (T, N) — > as 
N — ► oo. Well, Eq. [l^ cannot be valid at T — T c , because the lhs is finite while c(T c ) — oo (by hypothesis), i.e. the 
correction term has to be as large as the main term. This is not the only problem. Right at T = T c , N is the only 
parameter in hand, all others at their respective critical values. If we want -kC(T c , N) to diverge as N — > oo, we need 



C(T C ,N) 
N 



= C N p , (p>0),=>. C(T c ,N)~N 1+ p 



(13) 



which looks like a violation of Eq. (|5|) . We have assumed a power law form in the above equation , but many other 
possibilities exist. Our choice is motivated by the fact that a large number of systems (real or models) do show such 
power laws. 

Extensivity, ad nauseum, is a consequence of additivity: small pieces can be glued together to form a big piece with 
no change in property. In a sense boundaries can be ignored. This has to fail at the critical point. But how? 




FIG. 1. Schematic diagram showing the various paths in T vs N plane, around the critical point T = T c 



If we study the free energy, which cannot diverge, we might in principle write down an expansion of the type 
proposed in Eq. [12[ This expansion is in N for any T. From such a free-energy, we might compute the specific heat 
by taking derivatives and face again the problem of infinities on the right hand side at T = T c . Another way of saying 
the same thing is that the free energy at the critical point behaves as N 1+q with q < but the specific heat obtained 
by taking derivatives of the free energy with respect to temperature gets a different power of N. This means that 
N and temperature occur as a combination variable and not independently so that the limits N — > oo and T — > T c 
(taking derivatives) may not commute. The double limit needs to be considered carefully. 

The dilemma is shown schematically in Fig. 1. The expansion in N, Eq. (|l2]), is valid along path (a) ( T — > T c for a 
fixed N < oo) and path (b) ( thermodynamic limit for a fixed T), but the expansion is not valid in the neighbourhood 
of the critical point. The non-extensive behaviour gets reflected if the critical point (T = T Cl N^ — 0) is reached 
via a suitable path like (c) whose form is obtained in the next section. An isolated critical point does affect its 
neighbourhood as in any problem of non-uniform convergence. The neighbouring region (shaded in Fig [l]) is called 
critical region. 

Question remains: Can the specific heat be divergent? 



2 We have deliberately chosen an intensive variable. The case of the conjugate extensive variable is left as a problem. (See Ex. 
V.5.) 
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2. Comments 



• If the specific heat has the anomalous behaviour at 
T = T c , then the energy-density is also expected to 
have so. In general if the response function for a 
variable < $ > has a singular behaviour, the vari- 
able itself will also be singular. In the next section 
we see how these are related. 

• With Eq. (£f)in hand, we need to redo the argu- 
ment made in Se c II r egarding fluctuation. This is 
done in Section VIII. 



Question of divergence of specific heat: The answer 
is obviously "no" for a single particle system with 
a simple Hamiltonian (perfect gas, oscillator, two- 
level system ..). For a noninteracting collection of 
such simple systems, the free energy is always pro- 
portional to N and no criticality may occur even 
in the limit N — * oo. Interaction is needed and 
essential for criticality. 
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V. GENERALIZED SCALING 



The simple scaling of Eq. g needs to be modified to allow for the apparent non-extensive behaviour at a critical 
point. To show this we work in the ensemble with T as the variable (and all other parameters are kept fixed at their 
critical values). More variables will be taken up after that. We start with a generalization of scaling and show that it 
works. After working out a few consequences, we discuss the physical significance of the original scaling hypothesis. 

A. One variable: temperature 

Let us define t=(T — T c )/T c . Remember that this is an intensive variable and would not have required any scale 
factor on rescaling as per Eq. ^ or Ex. (0). In order not to clutter the equations with extra symbols, we consider a 
lattice problem for which N and V can be interchanged (e.g. the Ising model of Eq. (|3^)). Defining N = L d with L 
as a linear dimension of a (/-dimensional cube, we take a way out as 

F(t,L) = b~ d+Vi F(b~ Vt t,bL), (generalized homogeneity) , (14) 

from which Eq. (jE)) or Eq. (|J) can be recovered for yf = y t = 0. We first show that this generalized homogeneous 
function works and then, look at its consequences. Physical significance is discussed after that. 
Choosing b — L^ 1 , 

F(t, L) = L d L- y s T {jp^} ■ f(*, L) = b- y t f(b- y % bL), (15) 

where f(t,L) = jqF(t,L) is the free-energy per unit volume. We used T{x) — F(x,l). Finiteness of free-energy 
requires yf>0. 

For t ^ 0, demand extensivity. To get a linear dependence of F on N, we require^] T{x) ~ \x\ q such that 
~Vf + IVt — 0. This gives 

f(t) = km f(t,i) ~ \t\yf/y*, or f(t) = b-yff{wn). (ie) 

L — >oo 

Since the free energy per particle cannot be infinite even at T c , we need yf>0 and y t > 0. The last expression could 
also be written as f(t) = b~ y sl yt f(bt) by redefining b. 

If we are right at the critical point, the free energy density behaves as 

f(t = 0,L) ~ L- y f, (17) 

the "nonextensive" feature we are looking for. The generalized homogeneous function works. 

Exercise V.l Why cannot we choose F(T - T C ,L) = b- d+y fF(T - T c , bL) as a simple generalization? 
Hint: non-extensive everywhere. 

Exercise V.2 Is there any other generalization (other than Eq. ^ ) that could have been done? 

1. Comments 

A few things are to be noted in these manipulations. 



3 The sign ~ is used to denote the functional dependence on the variables on the rhs, coefficients etc may not be expficitly 
shown. Sign ~ is to be used to denote the feading term or terms with ail coefficients. 







A singularity of the free energy in the thermody- 



namic limit is manifest in Eq. 16 and the singular 
behaviour is a power law in the deviation from the 
critical point. These powers are called exponents. 

• We at the end recovered the much coveted exten- 
sivity, everywhere save the critical point, but at the 
cost of a scaling of the intensive variable. 

• What we are focusing here is the singular part of 
the free-energy. There could ( and will) be an ana- 
lytic piece that will not show any such anomalous 
scaling. 

• With a positive yt , we realize that the effective tem- 
perature is further away from the critical tempera- 
ture as the size is increased, if we want to keep the 
free energy per particle ( or per unit volume) the 
same. This desire to keep the same free energy per 
unit volume is consistent with the content of Eq. 



• A more dramatic result is that an argument 
L/jil^ 1 / 2 " can be thought of as a comparison of 



the length or size of the system with a character- 
istic scale £ of the system. This scale £ diverges 
when t 



as 



where 



1/Wt 



• A diverging length scale is the hallmark of a criti- 
cal point and in fact a critical point or criticality is 
defined as a point with a diverging length scale. 

• No, we have neither devised a way of bypassing 
the Legendre transformation problem nor ignored 
it. We chose the right ensemble (i.e. right thermo- 
dynamic potential) to do our analysis. It is con- 
ceptually helpful to use the ensemble that involves 
the intensive variables. For example, in first-order 
transitions (like, say, solid-liquid transition) there 
would be discontinuities in some extensive variables 
like the density, the energy density, the entropy 
( latent heat) etc, but the intensive variables re- 
main continuous. It is this fact that prompts one 
to use vertex functions in field theoretic analysis - 
but that's beyond the scope of these notes. 

• It still begs the question: What is this length £? 



2. Specific heat: power laws and exponents 

The specific heat is c - d 2 f/dt 2 . From Eq. ©, choosing b = \t\~\ it follows that c « /"(±1) \t\- 2+ nt/vt ^ w h ere 
prime denotes derivatives and ± stands for i^O. Note that the specific heat shows a power law singularity with an 
exponent 



y f — 

a = 2 -, defining a via c ~ \t\ a , or a 

Vt 



y lOgC 

hm — 

t-o log t 



(18) 



The exponent for divergence of c is same on both sides of the critical point at t = 0, though the amplitudes /"(±1) 
need not be the same. 

For L < oo, the specific heat behaviour is 



c(t,L) = L a l v C{tL x / v ), 



(19) 



making the L-dependence explicit as required by Eq. (|l3|). We leave it to the reader to establish the relationship 
between C of Eq. (|l|) and T of Eq. ( @. 



Exercise V.3 The two-dimensional ferromagnetic Ising model, Eq. in zero field (h = 0) has a 

logarithmic divergence of specific heat, c~ — ln|f|. How will the above formulation handle this? Note also 
that lnx = \im n ->o(x n — l)/n. 



3. On exponents: hyper-scaling 



Since free energy and energy-density cannot diverge, we need to have a < 1. This puts a limit 

Vf 1 
— > 1 or y f > yt or v > — . 

yt Vf 



(20) 



At the critical point ksTc is the energy scale, and therefore FJ (L d fcsT c ) dimensionally behaves like inverse volume. 
In absence of any other scale, one might expect 
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F(t = 0,L) 



Lr 



= d 



A similar argument for L 



k B T c N 

oo, and t ^ 0, would then give, with £ as the important length-scale, 
/ 



knT r 



t 



+du 



dv 



This relation between a and v involving d is called hyper-scaling. Putting all together 



kBT c 



= L- d T 



In general 



(21) 



(22) 



(23) 



(24) 



where yf could be different from d if some other length-scale plays an important role at the critical point. 

We did get the value of yj easily under certain conditions, but it is not possible to obtain v. Needless to say, 
with other variables, there will be more exponents. The values of these exponents are to be determined either from 
experiments or from statistical mechanical calculations with appropriate Hamiltonians. 

Exercise V.4 If there is no singularity for a finite system, then how can there be a diverging length scale 
in say Eq. (p3)? How to interpret this properly? Hint: asymptotic expansion (see Ex V.5 and Eq. 35). 



Exercise V.5 A problem from school algebra: Consider the sum S( x , N) = En=o forx<l. Show 

that for large N, and x close to 1, S(x, N) = NS(Nt) where t = 1 — x, and S(z) = (1 — e~ z )/z. Recover, 
from this asymptotic form, the exact result (i) for N — > oo for any x (Path b of Fig. 1), and (ii) for 
x — ► 1 for any finite N (Path a of Fig. 1). Compute numerically S(x,N) for various values of N and 
x < 1, plot S(x,N)/N vs Nt and compare with S(z) (Path c of Fig. 1). Note the strong deviation from 
the asymptotic scaling for small N or large deviation of x from 1 . This is actually a comparison of the 
partial sums of the infinite series around the singular point. 



4- Role of fluctuations: upper critical dimension: I 

From the form of the free-energy, the finite-size scaling of the energy at T = T c can be obtained, namely E/N ~ 
L -(i-a)/V_ Let 

us now go back to Eq. (^). The relative width of the distribution for E is given by 

_ ~ L a/2u L (1 ~ a)/u L~ d/2 = L (vi ~ d)/2 (2b) 
E ' v ; 

remembering that N = L d , and using Eq. ([m]). If now hyper-scaling is valid, yf — d, and see, the relative width does 
not vanish even in the limit of L — > oo. Therefore for such a case, fluctuation plays a crucial role. If, however, yf < d, 
then the relative width vanishes in the thermodynamic limit. We then get a situation where we have a critical point 
with a diverging length-scale and possibly diverging response functions, but the effect of fluctuations is not strong 
in the sense that the probability distribution remai ns sharp around the average value. Similar results are found for 



other extensive variables also as shown in Sec V C 4 



One thing becomes clear: the dimensionality of the system is very important and the lower the dimensionality 
the stronger is the effect of fluctuations. We might expect that for large enough d, this condition yf < d will 
be satisfied. One may then try to understand such critical systems by considering the average system ignoring 
fluctuations altogether. This class of theories is called mean-field theory. 

If there is a finite value of d — d u above which yf < d, then fluctuations can be ignored for all d > d u . This d u is 
called the upper critical dimension of the system. Hope would be that the effect of fluctuations can be studied in a 
controlled manner around this d u with e = d u — d as a small parameter. This is the basis of the e-expansion for many 
systems. 

We, however, caution the reader that the finite-size scaling in a mean-field theory ( or d > d u ) is more complicated. 
We are concentrating only on the fluctuation dominated case. 
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B. Solidarity with thermodynamics 



The problem with extensivity can now be understood in terms of a diverging length scale. In the strict sense, there 
is actually no violation of extensivity. The idea of adding small pieces to build a bigger one makes sense provided the 
smaller pieces themselves are representative of the bulk. For systems with short range interactions and small intrinsic 
or characteristic length scales, this is reasonable because the "small size" effects or boundary effects are perceptible 
only if the size is comparable to these lengths. These effects are small corrections that can be safely ignored. In case 
of a diverging length-scale or characteristic length-scales larger than the system size, smaller pieces cannot be added 
up. Thus at a critical point with diverging length-scales, the notion of adding up smaller pieces fails. The length-scale 
£ is the appropriate scale for comparison and so at a critical point the whole sample is to be treated as a single one. 

At any temperature t ^ 0, £ ~ \t\ ~ v is large but finite. If we take a block of size £ d as a unit, then extensivity 
requires that the free energy be proportional to the number of such blocks or blobs i.e. F = (N/£ d ) fo ~ N\t\ dv fo, 
with fo as the free-energy per blob. If fo is independent of £, then / ~ /o \t\ dlJ , recovering Eq. (p2|) . In case fo, the 
free energy of a blob of volume £ d , depends on £, then extra £ contribution is expected and one gets Eq. ( p4| ) with 
Vf + d. 

Exercise V.6 Anisotropic system: A particular two-dimensional system of size L x x L y shows the follow- 
ing scaling behaviour 

f(t, L x , L y ) » Xf (L x t v * , L y t"y) , (26) 

X being the size dependent prefactor. There are now two different length scales in the two directions. 
What are the length-scale exponents? Find the relation (hyper-scaling) among a,v x and v y . What would 

be the form of X if one takes L x x oo strips. (Ans: X — L"^" x .) What would be X for oo x L y strips? We 
assume that the strips do not show any phase transition. What would be the right combination variable if 
both L x and L y are to be used in the prefactor? (Ans: X = (L x 1 ^ v " + L y x l Vy )- a , Why?) 

Exercise V.7 Cross-over: Think of a cylindrical geometry. The system is finite in one direction (length 
L ) but infinite in the remaining d — 1 directions. There will be a critical behaviour in this geometry if 
d is not too small. Consider the d — 1 to d dimensional crossover as the finite length L is made larger. 
Consider various situations: (a) d < d u , (b) d > d u but d — 1 < d u , and (c) d — 1 > d u . 



C. More Variables: temperature and held 



Generalization of the scaling of Eq. ( |i4| ) to more variables is straight forward. Take a variable h which represents 
another intensive quantity like pressure, magnetic field etc, measured from its critical value. The critical point is now 
at t — 0, h — 0. The free energy can be written as 



N- 



F(t, h, L) = b- d+y t F(b~ v n, b- yh h, bL), "^>f(U h) = b y f f(b- yt t, b- yh h), 
so that the same series of manipulations done for h = would lead to 



f(T,h,L)= -F(t,h,L)=C Vf F± 



h L 

t 



L-l/u' L -y h 



(27) 

(28a) 
(28b) 



where T±(x, y) = F(±l, x, y). 

If tL Vt = const, then, as per Eqs. (28b) and (|2^), extensivity is nowhere to be found - that's curve (c) in Fig 1. 
The scaling form of Eq. ( |23| ) or ( 28b ) is known as finite-size scaling. 
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1. Comments 



• For a fixed iV, if t -> 0, then f - L"»/ . 

• For a fixed t, if N — > oo, then 

ft 



|t| 2 -«JT ± 



\t\Vh/Vt 



• We see the general feature: 

critical region "nonextensivity" 
size scaling Scaling. 



(29) 



finite 



• Once we take the diverging length-scale as the sole 
scale for the problem, the finite-size scaling form 
can be obtained from the bulk behaviour as well. 
A finite length L matters only when it is compa- 
rable to £. For L £, the above-mentioned blob 
picture is valid and for L < £ the whole system 
needs to be treated as a critical one. If the bulk be- 
haviour is like c(t) ~ \t\~ a , then the size-dependent 
cross-over is given by c(t,L) ~ \t\~ a C(L/^,h/^ Vh ) 



~ L a l v C(t/L~ 1 l v : h/L Vh ). (SeeEq. @.) Histor- 
ically, bulk scaling was proposed first and finite-size 
scaling came later on. 

• We cannot overemphasize the fact that h and t 
are completely independent variables. Yet in Eq. 

the free-energy in the thermodynamic limit de- 
pends on the single combination variable h/t y f' Vt , 
and not on h and t separately. Results of exper- 
iments (real or numerical) done at various values 
of h and t can be collapsed on to a single curve, 
unthinkable away from the critical region. 

• It seems there is one exponent too many for the 
bulk free energy density / in Eq. (p7|), yj/yt and 
Vh/yt would have sufficed. Yes, it would if we 
consider thermodynamics alone. But the path via 
finite-sized system showed us the necessity of an- 
other exponent, v, and so we keep all the three in 
Eq. A 



Exercise V.8 Take the Legendre transform of f in Eq. ( j£jj ) with respect to t. Discuss its scaling property. 
Caution: The scaling is not for the total entropy or entropy density. The special scaling is for the deviation 
from the critical value of the entropy at the critical point. This shows the difference of the scaling around 
the critical point and the simple scaling of thermodynamics. 



2. More variables mean more exponents 

The power laws we saw earlier tend to suggest similar behaviour for other physical quantities also. We define several 
such practically important exponents, but shall ultimately see that all of these can be expressed in terms of the three 
VfiVuVh introduced earlier. 

For concreteness, the magnetic language of the Ising model of Eq. (|h]) is used. Let us assume that there is a 
critical point for the Ising model (and there is one) so that the free energy is given by Eq, ( ^9j). We have already 
defined a as the specific heat exponent. We define (3 for magnetization, 7 for susceptibility, 8 for critical isotherm. 
Apart from these thermodynamic exponents, two other expoents are needed v ( already introduced) for length-scale 
and 77 for critical correlations. 

Remembering that the magnetization is the derivative of the free-energy with respect to h, one gets m w 
\t\(vf-Vh)/vt^(h/t Vh / v f). From this, see that for h = 0, the magnetization vanishes at the critical point as 

m ~ \tf, with f3=(y f - y h )/y t . (30) 

For a ferromagnetic transition, there is no magnetization in zero magnetic field in the high temperature phase (para- 
magnet), and therefore .F+'(0) = but .F-'(O) 7^ 0. A quantity like m that describes the "ordering" of the system 
is called an order parameter and /3 is the or der-pa rameter exponent. We repeat that phase transitions need not 



necessarily have an order parameter (see Sec. VII BP 



Susceptibility \ is the response of magnetization, and we see in zero field ( h = 0) (x — dm/dh) 

X~\t\~\ with J=(2y h -y f )/y t . (31) 

At the critical point (t = 0) in presence of a field ( critical isochore ), we get 

mRih 1/5 mo, 6 = y h /(y f - y h ). (32) 

Here the amplitude mo = lim a; _, 00 x~( y f~ Vh ^ Vh J-' ± (x). In general, 5^1. One expects a linear relationship ("linear 
response") between "cause" (h) and "effect" (m), but that turns out not to be the case at criticality. A linear relation 
can never give an infinite %! 
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A major consequence of the homogeneity of the free energy is the power law behaviour of various physical quantities 
(various derivative of free energy) and all of these are obtained from the three basic exponents. In case of hyper-scaling 
(i.e. yf = d) we have a further reduction and only two exponents are needed for a complete description of the critical 
behaviour of a system. 

For a thermodynamic description this looks enough, but we already saw the usefulness of a length-scale based 
analysis. We need to define the length-scale properly and in the process we will find a more useful critical exponent 

V- 



3. Comments 



• The blob picture can be used for susceptibility near 
the critical point. There are N/£ d blobs. Each blob 
of size £ d can be thought of as a critical object. The 
total susceptibility is given by N\ — {N/Z, d )x\Aah, 



where Xbiob is the total susceptibility of a blob. 
Finite-size scaling predicts Xbiob = £ d £ 7 / 1/ so that 
we obtain \ ~ \t\ " 7 

• Note that we require yh < yf < 2yh- 



4- Role of fluctuations: upper critical dimension:!! 



The role of fluctuations was analyzed earlier in the context of specific heat. Let us now reanalyze it for the 
probability distribution for magnetization M. The finite-size scaling behaviours of (total) magnetization and (total) 
susceptibility at T = T c are (M) ~ L d L~^l u and Nx ~ L d L~«l". The relative width AM/M of the probability 
distribution for M is £7/2^/3/^-^/2 = L ( Vf -d)/2 where Eqs _ Q and © have been used. 

It is reassuring that the conditio n for s harpness (or broadness) of the probability distribution for M is the same as 
obtained for specific heat (see Sec. V A4 ). The upper critical dimension is the same no matter which conjugate pair 
we use. 



D. On exponent relations 



Once the exponent identifications are made, with only two independent ones, it is possible to write down many 
relations involving the (in principle) experimentally measurable exponents. For example, by adding the exponents, 

a + 2/3 + 7 = 2. (33) 

If we take the free-energy density ~ hm and m = x _1 /i, then it follows that 2/3 + 7 = 2 — a, as in the above equation, 
without using any explicit formula. Another way of re-writing the above relation (—a = 2({3 — 1) — 7) suggests that 
the behaviour of specific heat (ch=o) is similar to (dm/dt^x^ 1 ■ In fact, thermodynamics gives us the formula 

C h~ C ™ = T (^j^j ( 34 ) 

where c x is the specific heat with x constant. Since specific heat is positive definite, it follows from Eq. (Q) for en 
with the intensive variable h held constant at the critical value h = that a + 2(3 + 7 > 2. 

Exercise V.9 Thermodynamic argument seems to indicate inequality rather than equality in Eq. (|5. 
Find out the conditions for which a strict inequality is expected. 
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VI. RELEVANCE, IRRELEVANCE AND UNIVERSALITY 



A few observations should not miss our attention. Since yt > 0, the combination variable Lt 1 ^* in, say, Eq. ( |l9| ) 



or ( 28b ) have different limits for t = and t ^ in the limit L — ► oo. This difference actually gave us the different 
L-dependent behaviour of the free-energy or the specific heat, or, as a matter of fact, any physical quantity we may 
calculate or observe. In the same way, if yh > 0, then in the bulk case, if t — > 0, i.e. as the critical point is approached, 
the combination or scaling variable h/t Vh / yt goes on increasing if ft, ^ 0, no matter how small it is. The behaviour is 
different if ft is strictly zero. Such a tendency of a parameter to grow also tells us that zero field and nonzero field 
behaviours are different. We call these variables relevant variables. 

In case yh < 0, then for t — > 0, the scaled variable is zero irrespective of its value. In such a case whether the system 
has ft, ^ to start with is immaterial. No wonder these are to be called irrelevant variables. 

There could also be redundant terms that don't matter at all, like a constant added to a Hamiltonian. We ignore 
them altogether. 

To be at a critical point, the relevant variables must be tuned properly to be at their critical values (e.g. t = h = 
in the magnetic example) because they take us away from the critical point. Irrelevant variables don't matter as such 
but they do play a significant role if we want to go beyond the leading behaviour. 

Special situations arise, if the scaled function shows a singularity as an irrelevant variable scales to zero. Such 
variables are then important for the critical behaviour, though they don't take the system away from criticality. Such 
a variable is called a dangerous irrelevant variable. 

Critical points are classified by the number of relevant variables required to describe them. An ordinary critical 
point requires two (for the magnetic problem temperature and magnetic field, for the liquid-gas transition temperature 
and pressure). If three relevant parameters are needed it is called a tricritical point and so on. 

Our analysis so far has been restricted to thermodynamic parameters like t, ft etc, but this can be extended to any 
parameter occurring in the Hamiltonian in a statistical mechanical approach. The starting Hamiltonian may have a 
large number of parameters based on the microscopic details of the system. But as we look at longer length-scales 
close to the critical point, all these parameters can be classified under the banner of relevance and irrelevance. By 
throwing away the irrelevant terms, for the leading behaviour, an enormous simplification ensues (e.g, all ordinary 
critical points will have only two relevant parameters, only difference may be in the numerical values of the exponents 
v and rj, etc.). It might then be expected that the numerical values of the exponents could be identified from certain 
basic symmetries etc of the Hamiltonian. This is the concept of Universality. A universality class would be described 
by the exponents and also the amplitude ratios of the various singular quantities on the both sides of the critical 
point. 

The idea of universality transcends the domain of critical phenomena. Whenever we are interested in properties 
on a scale much bigger than the underlying or microscopic scales, there seems to be a set of properties which are 
quantitatively same no matter what the microscopic details are. The exponents we have seen are just one such 
examples. Historically it was found that the shape of the coexistence curve near the liquid-vapour critical point is 
independent of the chemical composition and is also very similar to the critical behaviour of several magnets. Once 
the universality class of a system is identified, the set of universal quantities can be obtained by studying a simpler 
model system than the original one with all details. The simpler model is expected to focus on the relevant variables 
only or at most a few irrelevant ones. 

Universality is not just a set of exponents. In a scaling description, the amplitudes and even the scaling functions 
are independent of gross details except that the arguments of the functions may involve nonuniversal metric factors. 
In a finite size scaling there could be a dependence on the boundary conditions as well. 

Whether a variable is relevant or not is determined by its scaling exponent (e.g. yh in the previous case). In 



certain cases, one may get these by simple arguments (Gaussian model in Sec. IX) but in most cases these are to be 
determined. 

Exercise VI. 1 Can there be critical cases with only one relevant variable or no relevant variables at all? 
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VII. DIGRESSION: FIRST-ORDER TRANSITION AND TRANSITION WITH NO ORDERING 



A. A first-order transition: a=l 



What happens at the borderline of a = 1? The form of the free-energy tells us that the energy-density will have a 
discontinuity on the two sides of the singular point. Such a phase transition with a discontinuity in any first derivative 
of the free-energy is defined as a first-order transition. This value of a leads to v — 1/d. 

We see the possibility of a first order transition in the same framework developed for the critical point. But first- 
order transitions can be of other types also, and they may not necessarily have any diverging length-scales associated 
with it. One needs to be careful about it. 

As an example, let us consider a very simple configuration space for a magnet (a crude approximation to an Ising- 
type model). We replace the TV-spin configurations by two types only. A zero-energy ground state - all spins up or 
all down - two-fold degenerate, and an excited state of energy TVe with degeneracy 2 N . (This is a one-dimensional 
ferro-electric six vertex model, in disguise.) The partition function for this model is Z(x, TV) = 2 + (2a;) where 
x — exp(— fie). It is easy to see, by taking the TV — > oo limit, that there is first-order transition at x = 1/2 with total 
energy going from zero for x < 1/2 to TVe for x > 1/2. No problem with extensivity for x ^ 1/2. In the limit TV — ► oo, 
the specific heat is just a delta function at the transition point. For finite TV, the specific heat can be written in the 
form 

C(x,N) = N 2 In 2 (2 + e „ t)2 , (35) 

where t = 1 — 2x. For a (i-dimensional system, take TV = L d . What we now see from the scaling variable that there 
is a diverging length scale with exponent v = 1 j d. 



1. Comments 



• Fig 2 shows C(x,N) for various values of x and TV 
computed from the partition function and compared 
with the scaling function of Eq. |3^. This is an exam- 
ple of data collapse. 

• The thermodynamic limit for any given value of x 
(Path b of Fig. 1) comes from the tail of the scal- 
ing function of Fig 2. This plot amplifies the critical 
region of Fig. 1. 

• The TV 2 prefactor in Eq. ^ is consistent with an ex- 
ponent ctjv with a = 1. Remember that TV = L d . 

• The peak of the scaling function in Fig. 2 is not at 
the bulk transition point. 



Exercise VII. 1 What is this length scale 
in the context of this model? 



0.3 




N(l-2x) 

FIG. 2. Data collapse of C(x, N)/N 2 vs Nt for various TV. 
The solid line is the finite-size scaling function of Eq. B5l 



This simple model calculation can be generalized. One way of studying a first-order transition is to compute the 
free energies of the two phases independently ( like e.g. solid and liquid) and then finding the lower free energy curve. 
Taking these free energies per particle, /i and /2, as coming from a restricted sum over states of the full partition 
function, the total partition function can be approximated by Z — exp(— N (3fi) + exp(— TV/3/ 2 ), with the free energy 
per particle / = min(/i,/2) for TV — > oo. In general, /Vs are analytic function and the transition temperature is 
determined by fi{T c ) — f2(T c ), so that A/ = \fi(T) — f2{T)\ ~ \T — T c \ from a Taylor series expansion around 



1G 



T = T c . The partition function can therefore be written asZ = exp(-N0f) [1 + exp(-N/£ d )], where = c \Af\. 



This defines a length scale £ ^ \T — T c \ 1 l d . For first-order transitions, v = 1/d 



B. Example: Polymers : no "ordering" 

We give an example here of a phase transition for which there is no "ordering" unlike the magnetic case, and the 
interest in the critical point is because of finite-size behaviour. 

Polymers are long linear objects, abundantly occurring in nature. Let us take a single flexible polymer in a solvent 
( a single polymer of length N in d-dimensional space of infinite volume) . Depending on the nature of the solvent and 
the monomers constituting the polymer, there could be effective repulsion or attraction between the monomers. By 
changing temperature, one may go from a repulsion-dominated (monomer favouring solvent molecules) (self avoiding 
walk) to an attraction-dominated (monomer preferring monomer) phase. This is called a collapse transition. Such a 
transition for a single molecule can occur only in the thermodynamic limit of length N — > oo. But that's of no interest 
because the polymers are necessarily of finite lengths. The two phases in this example are described by the overall 
size of the polymer, e.g. by the mean square end-to-end distance (R 2 ) <~ N 2v . In a collapse transition the value of 
this size exponent v changes from v w 3/(d+ 2) to v = 1/d in rf-dimensions. This is an example of a phase transition 
or critical point which may not be viewed as ordering (unlike the magnetic system) - no less important though. The 
main features of the transition are reflected in the "finite-size" behaviour (with respect to TV) both at the transition 
and in the two phases, and the general analysis we have done so far is easily applicable here. 

Exercise VII. 2 What sort of a critical point is the collapse transition? Hint: How many relevant param- 
eters at a collapse transition? We need to make length going to infinity, we need to adjust temperature and 
if we take a solution of polymers, we need to make concentration equal to zero ( i.e. single chain limit). 
Ans: tricritical point. 

Exercise VII. 3 For a polymeric phase transition of the type mentioned above, show that the specific heat 
per unit length (or per monomer) can be written in a scaling form c = N 2 ^'~ 1 C{{T — T C )N^) 
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VIII. EXPONENTS AND CORRELATIONS 



A. Correlation function 

Consider the Ising model of Eq. ( [To| ) or ([ll]) and take M = J2i( s i)- The fluctuation-response theorem of Eq. (||) 
and translational invariance can be used to express the susceptibility (\ per particle) as 



N X = (k B T)-i - (s))(sj (*))> = N(k B T)~i J2 S 9(x) 



(36) 



s))) is the pair correlation function. The behaviour of the susceptibility is therefore 

bt 

($)2 



^ExS^W' where 



where = ((s - (s))(sa 

dependent on the behaviour of the correlation function. This relation is quite general as pointed out in section [II B 1 
By expressing the total quantity $ = 0(x), in terms of a local density, we have (<£> 2 
g H (x) = MO) - <0»(<Kx) - (0))). 

For simplicity, let's replace the sum in Eq. ( pq) by an integral so that x = fd d xg(x). For the Ising case, s is 
bounded and so is g. We also know that x diverges at least at T c . Only way this can happen is from the divergence 
of the integral or the sum - it is the large distance property of g(x) that controls the behaviour. We can conclude 
that, at least at T c , g[x) cannot be a short ranged function, but it has to decay to zero for infinitely large distances. 
In d dimensions, the decay at T = T c for large x is given by 



.-(d-2+fj) 



with r\ > (T = T C 



(37) 



For T T c , convergence requires that g(r) should decay sufficiently faster than this, better be of short-range character 
with a characteristic length scale £, e.g. exp(— x/£). The decay for any temperature can be written as 



1 



d-2+7) 



X 



(see comments below) 



(38) 



which goes over to the T — T c case if this length-scale £ diverges for T — > T c . The length-scale we saw popping out 
naturally for comparisons of the length of the system can now be identified as the correlation length, the scale for 
correlations. A precise definition^] of £ would be from the second moment of g(x) as 



J d d x x 2 g(x) 
J d d x g{x) 



\t[ 



(39) 



(For an isotropic system, the first moment vanishes by symmetry.) It is possible to define such scales via higher 
moments also. Under the assumption of one scale, all of these will have similar divergences at a critical point. 
Nevertheless, it is worth keeping in mind that there are cases where one may need to study higher order moments 
and different moments defining different length-scales. 

Exercise VIII. 1 Can 77 < 0? 



Exercise VIII. 2 Generalize Eq. 



for anisotropic cases like in Ex. V. 7. 



1. Comments 



The on-site term, i — j, with prefactor (ksT)^ 1 
in Eq. |3(j is the susceptibility of individual spins 
(or isolated degree of freedom), the Curie suscep- 
tibility. The correlation contribution is special to 
an interacting system. Note that the correlation 
vanishes for a noninteracting system. 



The fluctuation becomes long-ranged at the critical 
point. On the high temperature side (si) — ( no 
magnetization), and so the fluctuation correlation 
is the same as the spin-spin correlation. 

For low temperatures (T < T c ), the spin-spin cor- 
relation approaches a constant ((s) 2 7^ 0) for large 



There are in fact many ways of defining a length scale. We use a g(x)-based definition. 
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separations. This approach to the constant value 
is generally short-ranged, i.e., very rapid like expo- 
nential but becomes long-ranged at T c . 

• We get a long-range correlation even though the 
interactions are just short ranged ( e.g., nearest 
neighbour interaction in the chosen Ising model). 

• A finite size scaling for the correlation length itself 
would be £ ~ L. Quite often it is possible to write 
it as an equality with the amplitude determined 
by the known exponents. The amplitude surely 



depends on which correlation function is used to 
determine £, boundary conditions etc. 

• For historical reasons Eq. (|3^) is generally 
written in a slightly different form as g(x) — 

x d-2+r, D (§) exp(-x/f), with D(z) -> 1 for z -> 0, 
while D(z) ~ z (d-3+2»?)/2 for large z _ Such a 
form is required to make correspondence with the 
Ornstein-Zernicke theory of correlations, which we 
shall not discuss in these notes. 



B. Relations among the exponents 



Let's think of the pair correlation function. It decays rapidly once we are on a scale greater than the correlation 
length £. Close to T c , we can think of the system as blobs of highly correlated r egion s - the blobs are of size £ d in 
d dimensions. These are the blobs we introduced to salvage extensivity in Sec. |V B| . Inside a blob (r << £), the 
fluctuations are critical-like and at a simple level, a blob can be thought of as at T c , On a bigger length scale x » £, 
the blobs are independent. Basically, we are arguing that it is the correlation length that matters - all other length 
scales are unimportant. 

We use this simple picture for the susceptibility. We can cutoff the integral at x ~ £, and inside this region 
g(x) ~ x ~(d~2+v) _ xhe integral J dx x 1 ^ 71 ~ t; 2 ~ v . Using the temperature dependence of £, we get the temperature 
dependence of \ as I t \ ^ u ^ ri \ The net result is 

7=i/(2-»?). (40) 

This relation with the help of Eq. ( |3l| ) gives 

y h = (4i) 

This is a very important relation - it shows how the external magnetic field or the nonthermal relevant variable scales 
away from the critical point is determined by the decay of correlations at criticality. 
It is straightforward to see that a = v(2 — tje)- 



C. And it's correlations 



We now have the full identification of the three exponents: 



Vt = l/v, Vf = d, yh 



rf+2-r; 
2 



(42) 



All the nuances of critical phenomena are qualitatively and quantitatively expressed in terms of the correlation function 
and the exponents needed for it. It's a shift from a thermodynamic description to a purely statistical mechanical one. 
It is the correlations of the degrees of freedom that control completely the whole phenomenon. 

In this particular example of ordinary critical point, one of the exponents, r/ is really an exponent defined at 
criticality, while the other, v, is an off-critical one. An off-critical exponent helps in the description of the approach to 
the criticality. This can however be treaded for purely critical exponents. 



Exercise VIII.3 What ify f ^ d? 



1. Comments 
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Now that we have emphasized the role of corre- 
lation functions and expressed yh in terms of r\, 
why are we not doing the same thing for specific 
heat? Eq. (|]) tells us that specific heat is re- 
lated to energy-energy correlation function, and so 
shouldn't we define an exponent 77E? We could if 
we want to. But convince yourself, by using hyper- 



scaling, that yt = \/v 
use either v or rfc. 



{d + 2 - 77 E )/2. We may 



• One may generalize the analysis of Sec. VIII A for 
any local variable L({c/)(x)}) and the corresponding 
response function will be determined by the critical 
correlation of (L({(j)(xi)})L({(j)(x2)})} ■ 



In a correlation based approach, we need to compute the critical correlations of all possible combinations of the 
degrees of freedom (e.g. in the Ising case, spin-spin, energy-energy and so on). From these, the relevant variables can 
be identified. In the Ising criticality case, there are only two. In addition, for a thermodynamic problem, we need to 
know the finite size behaviour of the free energy at the critical point (or assume hyperscaling) . These three critical 
parameters then completely specify the leading behaviour even around the transition point. 

To repeat, a thermodynamic description would focus on yf,yt,yh for a critical point, while statistical mechanical 
description would focus on the set of purely critical exponents yj /yt, ?7e, and 77. The scaling relations we derived show 
their equivalence. 

Exercise VIII. 4 Can this happen: the spin-spin correlation function is given by Eq. (j?^) and the energy- 
energy correlation is given by gEE ~ x ~ tyd ~ 2+r,E ^ gE(x / £,e) where £e ~ \t\ ~ VE , but ve 7^ v? 



D. What's b anyway? : length-scale dependent parameters 

If we look at the free energy of Eq. ([l6]) or Eq. ( p7[ ) in the thermodynamic limit of L — > 00, a question may 
be asked, "what is b now?". Answer is in the correlation function. Under a rescaling x — > x/b, we see g{x,t) = 
b~( d ~ 2+r, > g(x/b, b Vt t), where the /j-dependence has been made explicit. This resembles the scaling of the free-energy. An 
interpretation of this equation is that if we change the scale of length measurements (actually the scale of microscopic 
details - though not apparent right now), the parameters of the problem and the concerned physical quantity get 
rescaled. Quantitatively, for x — > x' = x/b, we have t — ► t' = b yt t, h — ► h' = b Vh h, and /—>/' = b Vi f . 

For a given problem (i.e. t,h,.. fixed) such a scale transformation leads to a transformation of the parameters 
and the physical quantities. Repeated applications of this transformation would land us on a set of parameters 
which are functions of b. These are the scale dependent parameters - criticality is best understood in terms of these 
scale-dependent parameters. 

The scale-dependence of the parameters is best represented by infinitesimal scale transformation b — 1 + 51 (compare 
with Ex. III. 3). We may write down a differential equation for Eq. ( |l6| ) as 

(C(*) ^-Vt) / = °> ( 43 ) 

where Q(t) — dt/dl — y t t. This equation has the expected solution / ~ t y f^ Vt . 

This is a new way of looking at the problem. Instead of studying a system for various values of the parameters 
( "coupling constants"), we are studying it at various scales to see how it behaves in the long-scale limit, after all 
thermodynamic or macroscopic behaviour is for a large system. 

The equation derived just above, (Eq. (^), describes the flow of the free-energy as the scale is changed, and such 
equations are called flow equations. Any physical quantity will have a flow equation associated with it. However with 
only a few relevant parameters and a few independent exponents, not all flow equations are necessary. A fuller picture 
emerges from a renormalization group analysis. Actual flow equations turn out to be slightly more complicated and 
the simpler equation of Eq. (^) is obtained under special conditions of "fixed points" . The critical behaviour and 
the universality classes are ultimately linked to the fixed points. 

Exercise VIII. 5 Possibility of a dangerous fixed point?: A fixed point by definition is a point that does 
not change under rescaling. If a parameter say u has a stable fixed point at u — u* then u — u* can be 
taken as an irrelevant variable. Can this variable be a "dangerous irrelevant variable", or, in other words, 
can there be a dangerous fixed point? 
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IX. MODELS AS EXAMPLES: GAUSSIAN AND <j) 4 



At this point it is helpful to consider a few examples. The Ising model has been introduced in Ex. III. 6. Here 
we consider a continuum version of it. A naive continuum limit would give a Gaussian model for a field variable 

— oo < 6 < oo, 



H 



+ r)4> k ct>- k ] - h<j>k=o, (44) 



with a cut-off in real space, i.e. |x| > a or in momentum space \q\ < A ~ a , where a could be the lattice spacing. 
We shall assume a spherical Brillouin zone. A better continuum limit is the <j) A model, 



H 



Jd d x[^^) 2 + ^ 2 + u^}. (45) 



Exercise IX. 1 Starting from the Ising model on a hyper- cubic lattice of coordination number q = 2c?, get 
the continuum form given above. Show that r ~ T — T m with T m = qj . 

The reason for considering a continuum limit is to take advantage of dimensional analysis. Both sides of Eqs. 
and ^| being dimensionless and the right hand side expressed as a volume integral helps us in introducing a length 
based analysis in a natural way. 

That r = is a singularity is obvious from Eq. (B) because of the instability at r < 0. One can derive all the 
thermodynamic properties for the Gaussian model and convince oneself that there is a singularity at r = 0. We do 
not go into that. Here let us accept that r = corresponds to a critical point. 



In general, one would have a term of the type ic(V</>) 2 in Eqs. 44 and |4^. However, if c does not change sign, 
one may absorb it in the definition of <f> and redefine r and u. This has been done in those two equations. There 
are problems where c may change sign when external parameters are changed (Lifshitz point), and in such cases 
c needs to be kept explicitly. It would be necessary to keep c explicitly in case the field variable cannot be scaled 
arbitrarily, as e.g. if (j> is like an angular variable. 



A. Specific heat for the Gaussian model 

Using Gaussian integrations, one can compute the zero-field specific heat for the Gaussian model. The leading term 
is given by the integral 

2 (Stt)" J q (1 + q z Y 

where £ = r -1 / 2 , r — a 2 (T ~ T m ) and Kd is the surface integral of a d— dimensional unit sphere. The length-scale 
exponent is v — 1/2. 

Exercise IX. 2 Derive Eq. Show that t)e = d — 2. 

A dimensional analysis gives 

[cf>(x)} = i( 2 - d )/ 2 , [r] = L-\ [A] - L~\ and [u] = L d -\ [h] = L^ 2+d ^ 2 . (47) 

where L is a length scale. This simple dimensional analysis already identifies a diverging length-scale £ = r -1 / 2 , used 
in Eq. (^6|). Note also that the specific heat is a volume integral of the 4> 2 -(j> 2 correlation function in real space. A 
dimensionally correct form is therefore 

c = £ 4 - d C(u^- d , h^+^Z 2 , AO- (48) 

We take the Gaussian model first (u = 0) and no external field, h = 0. Then c = £ 4 ~ d C(A£) and Eq. [l6] is in this 
form. Take the limit £ — > oo. The behaviour of the specific heat now depends on C(z) as z — * oo. If this limit is 
finite, then the cut-off can be completely forgotten and £ plays the important role. In such a case the microscopic 
parameters are not important for the critical behaviour. This happens, as we see from the integral in Eq. (Eg|) for 
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d < 4, and a = (4 — d)/2, a value that satisfies hyper-scaling with v = 1/2. However for d > 4, the integral diverges 
at the upper limit and C(z) ~ z d_4 , leaving us with a non-divergent specific heat, i.e. a = 0. Hyper-scaling also gets 
violated. For h ^ 0, Gaussian integrations yield = (d+ 2)/2, consistent with dimensional analysis. The exponents 
are 



ijf = d for d < 4, but = 4 for d > 4 and v = 1/2 and 77 = 0, for all d . (49) 



In the above example we chose the specific heat because of its interesting behaviour for d < 4 and d > 4. Take 
Susceptibility. Dimensional analysis gives x ~ r . Convince yourself, by using Gaussian integrals, that this is so for 
all d. However the finite-size scaling behaviour will be different for d < 4 and d > 4. 



From the value of a or yf, (see Eq. Jl9j) and the sharpness criterion of Sec. VA4, we see that fluctuations can be 



ignored if d > 4. In a sense d — 4 turns out to be the upper-critical dimension for this model. 

Exercise IX. 3 Calculate the average energy of a Gaussian correlated blob. i.e. a blob of size £ d - do this 
by calculating the form of the energy for the Gaussian model and then putting r = 0. Study its behaviour 
for various d and then integrating once with respect to r estimate the free-energy /o of a blob (See Sec. 



V E ). Justify the violation of hyper-scaling for d > 4. 



B. Cut-off and anomalous dimensions 

The above warm-up exercise shows that the cutoff, a relic of the microscopic features of the model, cannot always 
be ignored even-though the relevant length-scales at which the phenomenon is taking place is much much larger than 
this. Dimensionality is also important. Dimensional analysis is not expected to yield any unique or useful relation if 
there are multiple scales in the problem. This is a very important point and we elaborate on this further. In case the 
cutoff can be ignored, the exponents are the same as predicted by dimensional analysis. This we see directly in the 
Gaussian model for d < 4 - the cut-off just didn't matter. 

Exercise IX. 4 Consider a finite Gaussian model with periodic boundary conditions in all directions. Ob- 
tain the finite-size scaling behaviour of the specific heat for various d. Note the discrete k sums with no 
k = mode. 

Exercise IX. 5 Do the finite-size scaling analysis for the zero-field susceptibility of the Gaussian model 
for various d. 

If we go back to the <^> 4 problem, then we need to treat Eq. ^|. Question now is what happens for £ — > 00. If C 
goes to a constant, this cutoff can be safely ignored. A general situation would be 

C(x,y,z) = z p C(xz Pl ,yz P2 ), as z -> 00. (50) 

Setting A = 1 , 

c _ £(4-d)+ P J( u £(<i-4)+Pi ) / l £y h +P2)_ ( 51 ) 
This has the expected scaling form when £ is replaced by In fact the same analysis can be done for £ itself, 

Z = r- 1 ' 2 X(ur^- d ^ 2 ,...,A/^), (52) 

and, in the r — > limit, the rhs of Eq. |5^ might (and would) pick up extra powers of r from the A-dependent 
argument]^]. This would then change the temperature dependent exponent of Eq. 5l| . 

Seen from a thermodynamic point of view, these extra powers of £ are remarkable because these seem to vitiate 
dimensional analysis. An additional scale is needed, and this comes from the small length-scale a ~ 1/A whenever 
there is fluctuations at all length scales. The ultimate exponent one observes are not what dimensional analysis 
based on £ as the important length-scale would have predicted, except for Gaussian-type models. (These dimensional- 
analysis-based exponents are called naive or engineering dimensions.) One needs to understand and explain the origin 



There will be a shift in the critical temperature also. We take r as the deviation from the actual critical temperature. 
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of the "extra" contribution like the p's, which are to be called anomalous dimensions. Once the role of cut-off is 
recognized, the discrepancy with dimensional analysis (which is infallible in any case) vanishes. 

Because of the long range nature of the correlations, any local fluctuation can affect regions away for it. This is the 
origin of scale-invariance or occurrence of power laws (scaling). However the occurrence of anomalous dimension is 
something more. In a Gaussian type model, the degrees of freedom can be decomposed into independent modes ( e.g. 
by going over to Fourier modes - "normal coordinates" ) and then we see the emergence of long range correlation in 
the long distance limit of q —> 0. The modes remain independent so that the fluctuation at one scale determined by q 
does not affect the fluctuations at other scales. For this reason dimensional analysis gives correct result. In contrast 
for a </> 4 type model, the modes for different g-values are no longer independent (coupled by the </> 4 term) and now a 
fluctuation at a short distance scale ( q close to the cut-off) can affect the fluctuations at longer scales even to q — > 0. 

The mode coupling allows a small fluctuation at some point in space to affect other points and the disturbance seen 
by any other point is the sum over all the paths that connect the two points in question. For low dimensions, this 
sum can have fluctuations and this leads to anomalous dimensions. In high enough dimensions, the availability of a 
large number of paths (phase space volume) helps in averaging out the effects. These two cases are separated by the 
upper critical dimension. 

• Note that u becomes irrelevant for d > 4 and it could be a dangerous irrelevant variable. The exponents in such 
cases would depend on the function also. 

C. Through correlations 

It is reasonable to expect that the scaling variable one sees in a given problem should be independent of the actual 
physical quantity one is looking at. Therefore in the limit of A — > oo, anomalous dimension like p2 for h should be 
the same for all quantities that depend on h. Let us take the example of the correlation function, in zero field, at 
criticality, 

9(x) = ^g(ux*- d ,Ax), (53) 

suppressing the arguments on the lhs. Now in the limit Ax — * oo, if g(...,z) = z~ v g(....), where the nature of the 
other arguments are not so crucial for us right now, we have 

9( x ) = ^2+^g(-) Mting A = I). (54) 

This r\ changes yt = (d + 2)/2 to = (d + 2 — rf)/2. Identification can therefore be made: p-2 = —r). 

Let us reanalyze the zero field critical correlation function g(r, A), where the A-dependence has been made explicit. 
Under a rescaling of all lengths x — > bx, we see g(x, t = 0, A) = 6~( d ~ 2 ) g{bx, 0, A/6). The scale factor for the correlation 
function picks out the exponent one would expect on dimensional analysis. However if a scale transformation is done 
that changes the longer length-scales but not the microscopic ones, i.e. £ — * 6£ but A — > A, then 

g(x, 0, A) = b-^+^gibx, 0, A) (55) 

In analogy with the naive dimension, we now define a scaling dimension which is the dimension one observes in the 
long scale limit keeping the microscopic lengths same. 

Denoting the scaling dimension of X by Sx from how X-X correlation function scales as in Eq. (]55|), while naive 
dimension by dx ( defined by dimensional analysis), we have <S^ = — (d — 2 + T})/2 while d$ = — (d — 2)/2. It is easy 
to check now that 

S h + 50 = -d = d h + d^. (56) 

This relation for the every pair of conjugate variables is very useful. The scaling behaviour of say Eq. ( pl| ) can 
then be interpreted as dimensional analysis but with scaling dimensions. Renormalization group transformation is a 
way of doing a transformation that scales the long lengths keeping the short ones same, thereby picking the scaling 
dimensions. 

Generalizing the above discussion, we may define the scaling dimension for any local quantity ?/>(x) from the 
critical auto-correlation function, i.e. the long distance decay of the ip-ij) correlation. The analogue of dimensional 
analysis would tell us that for the critical correlation of any combination of local variables ip±,ip2,..., we should have 
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(M0)M*2)-Mx n )) = \x2\- Sl - 8 '- 3n y (-,..., — ) . (57) 

where Si is the scaling dimension of ipi . 

Exercise IX. 6 Why is it that the dimensions add up to —d? When can it be something different? 

Exercise IX. 7 Use dimensional analysis to write free-energy /(r, u, h, A) and magnetization m(r, u, h, A) 
as 

f = C d K< 4 - d M [d+2)/ \K), (58a) 
and 

m = r (d - 2)/2 M« 4 ~ d , ^ (<i+2)/2 , AO, (58b) 
where f, M are unspecified functions. For A£ — * oo one expects anomalous dimensions: 

f = C d fH A ~ d+Pu ,M {d+2)/2+Ph ), (58c) 
and 

m = r (d-2)/2+PrnM( u ^-d+ Pu>h £(d+2)/2+ Ph y (58d) 



Express these new functions in terms of the functions in Eqs. fi58<\) and (58b). Using the relation m 
df/dh\ h=a , show thatp m = p h (= -rj). 
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X. CONCLUSION 



We attempted to give an introduction to critical phenomena especially the idea of scaling, its need and consequences. 
The emphasis is on diverging length-scales. In case of a diverging length scale, as at a critical point, the correlations 
become long ranged. Such a point is charact erized by (i) the decay exponents of the correlations ( which could have 
anomalous part like rj and t\e of Sec. VIII C), and (ii) the number and nature of the relevant variables at that point. 

Any problem that does not have any intrinsic or important length-scale would behave like a critical system. This 
absence of length-scales leads to power law decays of correlations and also power laws for other physical quantities. 
In low dimensions ( less than the upper-critical dimension, which could very well be infinite) fluctuations play an 
important role near or at the critical point and show up in the anomalous exponents. In such cases (i.e. d < d u ), the 
finite-size effects can be understood in terms of the finite-size scaling. 

Distinction needs to be made between scaling and anomalous dimension. Scaling is the rule for criticality, originating 
from a diverging length scale, while anomalous dimension is seen for fluctuation dominated cases. 

Renormalization group provides the proper framework for analyzing such phenomena. We feel it is worthwhile to 
motivate those ideas behind RG at the introductory level. Mean-field theory that ignores fluctuations could then be 
placed in the proper perspective as special (called dangerous) behaviour of certain irrelevant variables. 
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